New native QMD code in Geant4

Recently the number of hadron therapy facilities has been increasing rapidly. Some of them are designed for heavy ions and need a tool to simulate the passage of ions through matter. As well, the space engineering field requires good estimates of device damage caused by bombarding ions in cosmic rays. As a result, demands for detailed simulations of nucleus-nucleus interactions have increased. Quantum Molecular Dynamics (QMD) is the quantum extension of the classical molecular-dynamics model and is widely used to analyze various aspects of heavy ion reactions, especially for many-body processes, in particular the formation of complex fragments. A new C++ version of the QMD code has been developed for the Geant4 toolkit. Derivation of equations of motion and production of initial “ground state” nuclei are based on JQMD and R-JQMD. Recently proposed Lorentz covariant kinematics treatment in R-JQMD is also included. Final states of two-body collisions are calculated by the same scatter class which the Binary Cascade model uses. After the termination of the time evolution of QMD system, cluster identification carry out and identified clusters are regarded as a excited fragment nucleus and they are passed to excitation model in Geant4. Validation results represent a significant improvement over current Geant4 nucleus-nucleus collision models.